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Abstract 



It is well known in classical mechanics that, the frequencies of a periodic sys- 
tem can be obtained rather easily through the action variable, without completely 
solving the equation of motion. The equivalent quantum action variable appearing 
C I in the quantum Hamilton-Jacobi formalism, can, analogously provide the energy 

^ ' eigenvalues of a bound state problem, without having to solve the corresponding 

CT' Schrodinger equation explicitly. This elegant and useful method is elucidated here 

in the context of some known and not so well known solvable potentials. It is also 
^ , shown, how this method provides an understanding, as to why approximate quanti- 

^ I zation schemes such as ordinary and supersymmetric WKB, can give exact answers 

for certain potentials. 
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I. INTRODUCTION 

In classical mechanics, the Hamilton- Jacobi (H-J) theory is a well developed theory 
and provides an independent and often useful route for solving dynamical equations.^ 
In particular, for periodic motion, the action variable enables one to obtain the fre- 
quencies of a given system directly without having to solve the equations of motion 
completely. The quantum H-J theory has also been studied since the inception of 
quantum mechanics.^ It has been recently shown that, analogous to the classical pe- 
riodic systems, the quantum action variable can be profitably employed to arrive at 
the energy eigenvalues for potential problems, without obtaining the corresponding 
wave functions. ^'^ This is to be contrasted with the standard procedure to tackle 
bound state problems, where, the Schrodinger equation is solved both for the eigen- 
values and eigenfunctions. 

The whole approach of the quantum H-J theory to the potential problems is quite 
elegant, requiring only some knowledge of complex variables. Keeping in mind that 
a student with this background will be able to appreciate it, we have made this arti- 
cle quite pedagogical and self-contained. In Sec. II, we briefly outline the quantum 
H-J formalism and its connection with the Schrodinger equation and work out the 
familiar harmonic oscillator example exphcitly. Section HI is devoted to the study 
of some known and not so well known potential problems to elucidate the power of 
this method. In Sec. IV, we show why approximate quantization schemes such as, 
WKB and supersymmetric (SUSY) WKB, give exact answers for certain potentials. 
We end this paper with some discussions and concluding remarks. For convenience, 
a table giving relevant information about the potential problems has also been in- 
cluded. 
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II. QUANTUM HAMILTON- J AC OBI FORMALISM 

Quantum H-J formalism has been developed along the lines of the classical H-J 

theory since the beginning of quantum mechanics, by the pioneers of the field. ^ 
In fact, this approach was christened as the "Royal road to quantization" in the 
early days of quantum mechanics.^ In 1983, Leacock and Padgett^'^ demonstrated 
that the quantum H-J formalism can yield the exact eigenvalues for the potential 
problems provided boundary conditions are applied judiciously. In what follows, 
after elaborating on the connection of this method with the standard text book 
approach to bound state problems, we solve the harmonic oscillator potential, as an 
illustration. 

In the conventional approach to non-relativistic stationary state problems, one solves 
the Schrodinger equation 

= (l^ + V{x)^ ij = Ei: , (1) 

for the eigenvalues and eigenf unctions. In the quantum H-J formalism, the postu- 
lated quantum H-J equation, 

hd'^W(x,E) fdW(x,E)Y hdp(x,E) ,^ 

= A(x,E) , (2) 

replaces the Schrodinger equation as the dynamical equation. 
Here, W{x, E) is the quantum characteristic function and 

Pi^, E) = ' , (3) 

is the quantum momentum function (QMF) and Pc{x, E) is defined to be the classical 
momentum function: 

. (4) 
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This can be thought of as a manifestation of the correspondence principle or as a 
boundary condition on the QMF. As will be seen explicitly later, the above condition 
helps in determining p{x,E) unambiguously. The quantum characteristic function 
W{x, E) is related to the energy eigenvectors in the coordinate representation as, 

iW(x,E) 

ip[x,E) =< x\E >= e R , (5) 

and in the same representation, 

.9 ,^ dW 
< x\p\E >= -ih— < x\E > = — — < x\E > 
ox ox 

= p{x,E) <x\E> . (6) 



Thus one gets. 



hi dip{x, E) 



It is straightforward to check that the Schrodinger equation goes over to the corre- 
sponding quantum H-J equation under the above substitution and vice versa. 
The quantum analog of the classical action variable is defined as 

J (E) = {l/27r) (f dxp{x,E) . (8) 
J c 

Here, C is a counter clockwise contour in the complex x— plane, enclosing the real 
line between the classical turning points. The turning points between which the 
classical motion takes place, are the real values of x, for which pl{x,E) vanishes. 
The wave function is known to have nodes between the classical turning points. 
These nodes correspond to poles of the quantum momentum function. To see this 
clearly, near a zero of the wave function, located at Xq, we write, 

^{x - xo)0(x) . (9) 

This imphes, 

p^- -^ + ... . (10) 

I X — Xq 
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It is thus seen that, p{x, E) has a first order pole at xq with residue —ih. One can 
also verify the correctness of Eq. ( p!OD directly from the quantum H-J equation. 
Substituting Eq. (|10D in Eq. (D, one sees that the contributions of these poles from 
E) and —ihdp{x, E)/dx cancel each other only if these poles are of first order, 
each having the residue —ih. The first order poles are of quantum mechanical origin 
and their positions are energy dependent, being the same as the zeros of the energy 
eigenf unction. Just as the zeros of the wave function change their positions with 
energy, so do the location of the corresponding poles in the QMF. These poles will 
be referred to as the moving poles. The quantum H-J equation shows that, p{x, E) 
can have singularities in the complex x— plane, other than the moving poles on the 
real axis. These singular points correspond to the singular points of the potential 
term V{x). The locations of these singularities are, quite obviously, independent of 
energy; these poles will be called fixed poles. 

For a given energy level the quantum number '?7,' equals the number of nodes of the 
wave function and hence it counts the number of moving poles of p{x, E) inside the 
contour C appearing in the definition of the quantum action variable given by Eq. 
(H). The residue of each of these poles is —ih. Hence, we have 

J{E) =nh , (11) 

as an exact quantization condition. 

Even though this approach appears similar to that of the familiar WKB scheme, 
it is worth pointing out that, Eq. (^^, when inverted for E reproduces the exact 
quantized energy eigenvalues. Although a priori, the location and the number of the 
moving poles are not known, a suitable deformation of the contour in the complex 
plane and change of variables allows one to compute J{E), for many potentials, in 
terms of the fixed poles whose locations and residues are known. In what follows, 
explicit examples are given to clarify this method. 
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Examples: 

1. Harmonic oscillator 

The quantum H-J equation for the harmonic oscillator problem with V{x) — 
is, 

9 hdp(x,E) ,^ muj^x^ . 9 , „, 



The turning points, determined from pl{x, E) = 0, are —Xi = X2 = +^2E/ (ma;^). 
The quantization condition is given by 

J(E) = (l/27r) (f p(x, E)dx = nh . (13) 

Here C is the contour enclosing the moving poles between the two turning points xi 
and X2 (see Fig. 1). Noticing that, there is only one fixed pole of p{x, E) dX x ^ oo, 
to evaluate J{E), one considers an integral /r^ over a circular contour having 
radius R and oriented in the anti-clockwise direction. The QMF has no singular 
points between Vr and C. Hence, for this case, J{E) coincides with Iyji'. 

h.^m . (14) 

For the evaluation of the contour integral Itj^ , one makes a change of variable x — 1/y 

to get, 

Ir^ = (l/27r) / dxp{x, E) 

= (l/27r) / dyp{y,E)/y^ . (15) 



Here, p{y,E) — p{l/y,E) and the counter clockwise contour 70 encloses only one 
singular point in the y plane, i.e., the pole at y = 0. The corresponding contour 
integral can be straightforwardly calculated. Note that there is no negative sign 
before the integral; the direction of the contour changes sense under this mapping, 
which is compensated by the negative sign coming from the integration measure. In 
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this example J{E) and Jr^ are equal, though the relation between J{E) and Jr^ 
will change from one example to another, the method of computing Jp^ remains the 
same for all examples. 

The quantum H-J equation written in the y variable becomes 

p[y,E) + ihy — — — = 2m(E - ^p-) = . (16) 

To calculate the contribution of the pole at ?/ = 0, p{y, E) is expanded in a Laurent 
series as, 



oo ^ h 

p{y,E) = Y.any'' + Y.7-, ■ (17) 

n=0 (3=1 y 

Making use of the above expansion of p{y, E) in Eq. (]T5D, one notices that, the only 



non-vanishing contribution comes from the coefficient ai of the linear term in y. 
In the next step, substituting E) in the quantum H-J equation and comparing 
the left and right hand sides, it is found that, 6g = for g > 1. On equating the 
coefficients of the different powers of y, we have, 

hi = -m^uj'^ , (18) 
2ao6i = , (19) 
-ihhi + 2ai6i + = 2mE . (20) 

From Eq. (|l^), one finds hi = ^imu. This ambiguity in sign for hi can be removed, 
if we apply the boundary condition given by Eq. (^). In the convention followed 
here, the classical momentum function is defined such that Pc{x, E) = +i\pc\ on the 
positive real axis."^ In the limit ?/ — > 0, (x — > oo), pc ~ imuj/y and therefore from Eq. 
(|T7p it follows that hi = imu. From Eq. (|T^, we then have Oq = 0. Substituting 
the value of hi in Eq. (|20D one gets oi = {2E — fiuj) / {2iuj). Plugging Eq. ([TtD in Eq. 
( p!5|) and noting Eq. ([l^) gives. 



J{E) = Jr, = ^al = . (21) 
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Thus the quantization condition J{E) = nfi when inverted for gives 

E = (^+]^nw . (22) 
2. Harmonic oscillator on half line 

The quantum H- J equation for this case is obviously the same as that of the previous 
example, apart from the condition that the potential is oo at x = 0. This forces us to 
assume a fixed pole for p(a;, E) in the complex x plane at a; = and this also serves 
as one of the turning points. The other turning point is located at ^2 = \jlE l{mw^\ 
We have, 

J{E) = (l/27r) I p{x)dx , (23) 
J c 



where C is the contour enclosing the moving poles between zero and 2E / {mui'^) . 
To evaluate J{E), we consider a contour integral /p^ over a circle of radius R in 
the complex plane; R is taken large enough to enclose all the singular points of the 
QMF inside it. This contour integral, 

/r« = (l/27r) / p{x)dx , (24) 

gets contributions from the singular points of p{x, E) inside Vr. These are, (i) fixed 



pole at a; = 0, (ii) moving poles between and X2 = 2E / {mui'^) and (iii) moving 
poles on the negative real axis between and x^ — —^2E/{muj'^). The poles on the 
negative real axis arise due to the symmetry x — > —x of the problem. Hence, 

7r, = J{E) + Ic, + A, , (25) 

where I^^ is the contour integral for the contour 71 that takes into account the 
additional pole at a; = and Jp^ is the contour integral for the contour F/j enclosing 
all the moving and fixed poles of p(a;, E). is the contour integral for the counter 
clockwise contour Ci, enclosing the moving poles on the negative side of the real 
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axis (see Fig. 2). However, under x —x, the turning points xi and X2 and the 
moving poles get interchanged. Thus, 

Ic, = J{E) . (26) 

Substituting the above in Eq. (|^) we have, 

/r, = 2J{E) + J,, . (27) 

Observing that, it is only the hi/x term in the Laurent expansion of p{x,E) that 
contributes to the contour integral /^^, the relevant relation involving bi following 
from the quantum H-J equation is, 

bl-^h = . (28) 

From the two solutions hi = 0; —ih, the former is discarded, since it does not give 
rise to the required singularity at x = 0. Hence, using the second value, 

I^,=h . (29) 

/r^ is evaluated as in the previous section, and is found to be exactly the same. 
We now substitute the values of /^^ and Jr^ from Eqs. (^) and ( pi]) in Eq. (pT]). 
The quantization rule J{E) = nh, when inverted for E, gives the energy eigenvalue 
equation; 

E = ({2n +l) + -)huj . (30) 



2. 

III. OTHER SOLVABLE EXAMPLES 

In this section, the quantum H-J method will be applied to a set of trigonometric 
and hyperbolic potentials.^ These potentials have attracted considerable attention 
in the recent literature due to the fact that they can be solved algebraically using 
the techniques of SUSY quantum mechanics.^ As will be seen explicitly later, SUSY 
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shows an alternate way of selecting the correct solutions of the quantum H-J equa- 
tion. Suitable exponential mappings are used in the analyses of the trigonometric 
and hyperbolic potentials. Here, the example of the Eckart potential will be worked 
out in detail and the relevant information about the rest of the potentials will be 
given in Table I. We also work out the infinite square well potential in this section 
which has been deliberately tackled at the end because of its nontrivial nature in 
context of the present formalism. 
1. Eckart potential 

The Eckart potential, by suitably adjusting the ground state energy, can be written 
in the form, 

V{x) = uj''{x) 



2/^\ ^ duj{x) 



Phn dx 
where 

uj{x) = -Acoiliax + B/A , (31) 

is known as the superpotential in the literature.^ Given the ground state wave func- 
tion ipQ, in SUSY quantum mechanics, 

y/2mtpo ox 

The QMF, p{x,E), becomes equal to i\/2muj{x) for E = 0. The quantum H-J 
equation for the Eckart potential is given as, 

p^{x,E) + - ^^^^'^^ =2m\E - A^ - ^ - A{A - -^)cosech^ax + 25 coth ax ) 



i dx \ A^ \j2m' 

(33) 

where x lies on the half line. To simplify the analysis, we use the mapping 
y = exp(Q;x) and the corresponding equations (PT|) and (RSF) become 



y — 1 A 
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and 



9?/ \ {y^-iy (y'-l) 

(35) 

respectively. Here p{y, E) = p(\ogy/a, E). The right hand side of Eq. ( P5| ) equated 
to zero has four solutions, out of which two turning points are in the physical regime. 
The quantization condition in the y variable takes the form, 

J{E) = %{y, E)=nh . (36) 

zvra Jci y 

It is clear from Eqs. (|35|) and (|36|) , that the integrand has singularities at y = 0, ±1. 



As before, we shall now consider a circle of radius R, which is such that all the 
singular points of the integrand in Eq. (^) are enclosed. Then 

Ifr = Ici + Ic2 + Ai + + -^73 • (37) 

where = J{E) and are the integrals along the counter clockwise contours Ci 
and C2 enclosing the classical and non-classical turning points respectively. /^^, I^^ 
and are the integrals along contours 71, 72 and 73 which encloses the singular 
points at y = l,y = —1 and y = respectively (see Fig. 3). It may be noticed 
that, the symmetry y — > —y in Eq. ( p^ interchanges the turning points in the 
non-classical region with those in the classical region. Thus, 



(38) 



Therefore, from Eq. (|37|), we have 

2 J(E) = Jr, /,3 . (39) 
To find the contribution for the pole at y = 1, one expands 

p{y,E) = j;^^^ + ao + a^iy -!) + ■■■ , (40) 
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and substitutes the same into Eq. (^). Comparing the coefficients of — 1)^ on 
both sides gives, 

—iah ± i(ah — 2\/2mA) , 

&i = ■ (41) 

Similarly for the pole at ?/ = — 1 an expansion in powers of + 1) is sought in 
E) and we have 

. icxh ± iiah — 2\/2mA) , 

&i = ^ ^ , (42) 

where 6'^ is the coefficient of the + 1) term in the above expansion. To find the 
residue of the integrand at the pole at ?/ = 0, we expand i?) as 

E) = ao + aiy H , (43) 

and comparing the coefficient of the constant term gives 

al = 2m[E - - — -2B) . (44) 

For the calculation of Jr^, one more change of variable in the form of y = Xjz is 
sought, so that the singularity at y oo is mapped to the one at 2 = 0. Proceeding 
as before, the coefficient of the constant term in the expansion of j)(z, E) in powers 
of z is, 

o!l = 2m{E -A^- — + 2B) , (45) 

where p{z, E) = p{l/z,E). To find the spectrum for the Eckart potential, one has to 
choose correct signs for the above coefficients, by appropriate boundary conditions. 
The procedure as originally suggested by Leacock and Padgett is a bit complicated 
to implement. Instead, we adopt an alternate method, where a similar expansion of 
the superpotential uj{x) at the given poles is compared with the E limit of the 
above coefficients, there by fixing the correct sign. This is because Eq. (^) implies 
that p{x,E) = i\/2muj{x) at zero energy. Therefore from Eqs. (^T]), (|42|), (|4^) and 
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51), we have, 



al^, = V2^A , (46) 
aJ^2 = V2mA , (47) 
al^, = -id2m{E - A'^ - — -2B) , (48) 



air, = ^2miE-A^-^ + 2B) . (49) 

Taking into account the contribution of the two identical contours enclosing the 
moving poles in the complex ?/— plane, we find 

Ir, - E ^> = 2ri^ . (50) 

p=i 



Solving for E: 



52 52 , nah 

nah 



^ V 2m ' 



The calculation of eigenvalues for other SUSY potentials proceeds along the same 
lines and the results are summarised in Table I. 
2. The square well potential 

Square well potential of width L is the simplest example of one-dimensional motion 
where a particle of mass m experiences a potential, 

V{x) =0 for < X < L , 

= 00 for X < and x > L . (52) 

The nodes of the eigenfunctions and hence the moving poles of p{x, E) are located 
between x = and L. Let C be a rectangular contour enclosing all the moving poles 
(see Fig. 4a). Then the quantization condition is given by 

J{E) = ^j>p{x,E)dx = nn , (53) 
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where the QMF obeys the quantum H-J equation, 

p\x,E)-ih^^!^^ = 2mE . (54) 

We now use a mapping z = exp((27rix)/L); the contour C in the x— plane (Fig. 4a) 
is mapped into the contour PAQRBSP in the 2— plane . The moving poles get 
mapped on the middle arc of unit radius (Fig. 4b). The quantization condition is 
now given as 

where Ci is the contour PAQRBSP of Fig. 4b. The quantum H-J equation written 

in terms of the new variable z is, 

_n, 2nhz dp(z, E) ^ , , 

p\z.E) + — ' =2mE , 56 

L oz 

where p{z, E) = p{L \og{z) / {2ni) , E). The boundary condition that, the wave func- 
tion vanishes at x = and x = L gives rise to a pole in p{z, E) at z = 1. Let 7 and 
F be the inner and outer full circles of radii OA and OB respectively, both taken 
in the anti-clockwise direction. The integral in Eq. ( ^5]) can be written in terms of 
integrals over F, 7 and the contour P'S'R'Q'P' enclosing the pole at z = 1. Thus, 
we get 

j^^^^J^/rp{z^_rp{z^_r p{z,E)dz\ ^^^^ 

Ain"^ \Jr z z Jp's'R'Q'p' z J 

The first integral is computed by changing variables from z to y = 1/z, as was done 

earlier for T^. The last two integrals in the above expressions are calculated as usual 

by doing the Laurent expansion of p{z, E) in powers of z and z — 1 respectively and 

we have, 

/^^^^ = -2.^V2^ (58) 
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rpiz^ = 2mV2^ (59) 
p{z,E)dz in^ih 



IP'S'R'Q'P' z L 
Substituting the above in Eq. (0), and solving for E we get, 



(60) 



^=(l^)<"+^>' • " = ■ (61) 

IV. EXACTNESS OF ORDINARY AND SUSY WKB APPROXIMA- 
TION SCHEMES 

Semiclassical quantization schemes hke WKB have been very useful since the early 
days of quantum mechanics, in getting the approximate spectra of bound state 
problems. Interestingly, for certain potentials, these approximation schemes give 
exact answers. It is well known that WKB quantization condition, 

rx2 I 1 

/ J2m{E - V{x))dx = {n + -)-Kn , (62) 

reproduces the exact eigenvalues for the harmonic oscillator potential. Here, Xi and 
X2 are the two classical turning points with xi < X2 and n takes positive integral 
values. 

With the advent of SUSY quantum mechanics it was found that, for a potential 
V{x) = uj'^{x) — :^7^^^^~5 for which the ground state energy is zero; a variation of 
the WKB approximation, 

y""' y^2m(E - uj^x)) = nhn; n = 0,1,2, ... , (63) 

gives exact eigenvalues for a large class of potentials. Here xi and X2 are solutions 
of E — uj"^ = with Xl < X2. This is the well known SUSY WKB approximation 
scheme. Since quantum H-J formalism gives exact results, and is very similar to 
these schemes, it can be used to gain an understanding as to why SUSY WKB 
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approximation reproduces the exact answers for these potentials. 
For the harmonic oscillator, where WKB is exact, it can be shown that WKB and 
SUSY WKB approximations are identical. Therefore, we concentrate only on the 
SUSY WKB scheme here. In order to compare SUSY WKB with the quantum 
H-J quantization condition, the line integral will be converted to a suitable contour 
integration over the complex x— plane, with a counter clockwise contour C, enclosing 
the turning points Xi and X2 .^^ Writing, 



the contour integral in Eq. ( |64D can be evaluated, for all the potentials tabulated in 



Table I. Amazingly, for all these cases, the singularity structure of \j2m{E — tu^), 
other than the branch cuts, matches exactly with that of the fixed poles of E) 
in the quantum H-J formalism; the location of the fixed poles and the corresponding 
residues are identical for both the cases. As has been seen earlier, these poles and 
their residues determine the eigenspectra completely. Hence, it is not surprising that 
SUSY WKB gives exact answers for these potentials. 
V. CONCLUSION: 

We have explicitly worked out the eigenvalues for several solvable potentials in one 
dimension, using the quantum H-J method. Table I contains all the relevant details 
about the steps involved in the calculation. Interestingly, the fixed poles of the 
quantum momentum function, which are not of quantum mechanical origin, in con- 
junction with the boundary conditions, completely determine the spectra for these 
examples. The main effort involved in use of this scheme, lies in selecting the correct 
roots for the residues needed. This was done, first, by using the boundary condition 
on the QMF E), viz., p — in the limit — ^ 0. It was then shown that for 
the SUSY potentials, comparison of the answers obtained from the quantum H-J 
for E = {] with that obtained from the superpotential also reproduced the correct 




(64) 
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roots. 

For the potentials, where certain approximate quantization schemes are exact, this 
approach has provided some useful insight: It is the similarity of the singularity 
structure of p{x, E) and ^2m{E — o;^) in the non-classical regions of the x— plane, 
namely the matching of the poles and the residues that is responsible for this exact- 
ness. 

Although, we have not dealt with them here, for the non-exact potentials, it is the 
inability to deform the contour appropriately because of the presence of other poles 
and branch cuts in the complex x— plane that prevents an exact solution in this 
approach. Since, exactly solvable potentials are few and far between, this method 
can be potentially useful to construct new ones. 
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Table I : Hyperbolic and trigonometric potentials. 

The mapping used for hyperbolic potentials is y = exp(Q;a;), while 

for the trigonometric ones y — exp{iax) is being used; a is real and positive. 



Name of 


Potential 


Fixed al-yj, 


Eigenvalue 


potential 




poles at 





Scarf II 



Aah \ 
ah 



X 



(hyperbohc) sech^ ax + B{2A + ^) x 
sech ax tanh ax 





i 

—i 
oo 



V2^{iB - A) 
-V^{iB + A) 
i^2m(E - A^) 



A' -{A- 



2m' 



Rosen - 
Morse 11 



+ ^2 

-A{A + 



ah 
\/2m 



sech^ ax 



(hyperbolic) +25tanh ax 





i 

—i 
oo 



-ipm{E -A^-^ + 2B) A^ + B^A^ 
-V2^A -{A - ""^^2 

-\/2mA 

i^2m{E -A^-^-2B) 



B'^HA-^f 



Generahsed 

Poschl- 

Teller 



2m 



A^ + + A^ + 
cosech^ ax — B{2A 
coth ax cosech ax 



IX 

ah 
\/2m 



(hyperbohc) {x > 0) 




1 

-1 

oo 



-i^2m{E - A^) 
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Figure Captions 

1. Figure 1. Contour for the harmonic oscillator problem. 

2. Figure 2. Contour for the harmonic oscillator problem on half line. 

3. Figure 3. Contour for the Eckart potential problem, using the mapping 

y — exp(Q;a;). 

4. Figure 4a. Contour for the square well problem, in the x— plane. 

5. Figure 4b. Contour for the square well problem, in the plane. 



